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Abstract 

. It has been known for some years that entanglement entropy obtained from partial trace does not 

00 

', provide the correct entanglement measure when applied to systems of identical particles. Several 

' criteria have been proposed that have the drawback of being different according to whether one 

■ is dealing with fermions, bosons or distinguishable particles. In this Letter, we give a precise and 

m 

mathematically natural answer to this problem. Our approach is based on the use of the more 
general idea of restriction of states to subalgebras. It leads to a novel approach to entanglement, 



^ ' suitable to be used in general quantum systems and specially in systems of identical particles. This 

settles some recent controversy regarding entanglement for identical particles. The prospects for 
applications of our criteria are wide-ranging, from spin chains in condensed matter to entropy of 
black holes. 
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INTRODUCTION 



The study of subsystems of a quantum system is of paramount importance in many 
branches of physics. In quantum information it enters in the analysis of local operations 
performed by different parties of a multipartite system. In statistical physics it enters in 
the very definition of the different ensembles since this involves considering a given physical 
system as embedded in a bigger one. In the physics of black holes, the distinction between 
accessible and inaccessible regions of space-time plays a crucial role for the study of black 
hole entropy. Indeed as pointed out in [1], the coupling from outside to inside the horizon 
is very strong, while the reverse coupling is nonexistent! In all these situations partial trace 
is the preferred tool to extract physical properties of the given subsystems. Nevertheless, it 
is well-known that in some cases of great physical interest like systems of identical particles 
the use of partial trace leads to contradictory results. 

In this Letter we provide a resolution of such contradictions which turns out to be of 
general application. We show that, by treating observables and states on an equal footing, 
a generalized notion of entanglement emerges. A relevant consequence is that the entangle- 
ment measure that naturally arises in this algebraic approach is shown to be easily computed. 
Our approach thus opens up a wide range of applications, from condensed matter systems, 
like spin chains and anyonic models, to black hole physics. 

For bipartite systems contradictory results due to partial trace are explicitly seen to 
appear in the computation of entanglement measure for identical particles systems. In spite 
of the numerous efforts to achieve a satisfactory understanding of entanglement for systems 
of identical particles, there is no general agreement on the appropriate generalization of 
concepts valid for non-identical constituents IJIsl. That is because many concepts are 
usually only discussed in the context of quantum systems for which the Hilbert space H is a 
simple tensor product with no additional structure like, for example, 1-i = Ha ^'Hb- In this 
case, the partial trace pA = Tr^ for \tp) G "H to obtain the reduced density matrix 

has a good physical meaning: it corresponds to observations only on the subsystem A. 

In contrast, the Hilbert space of a system of identical bosons (fermions) is given by the 
symmetric (antisymmetric) A^-fold tensor product of the single-particle spaces. The conse- 
quence is that any multi-particle state contains intrinsic correlations between subsystems 
due to quantum indistinguishability. This, in turn, forces a departure from the straightfor- 



ward application of entanglement-related concepts like singular value decomposition (SVD), 
Schmidt rank or entanglement entropy. 

We propose here an approach to the study of entanglement where the notion of partial 
trace is replaced by the more appropriate notion of restriction of a state to a subalgebra^\. 
This approach is based on the well established GNS construction[10]. It allows us to mean- 
ingfully treat entanglement of identical and non-identical particles on an equal footing, 
without the need to resort to different criteria according to the case under study. 

The usefulness of our approach wll be d.splaved i,r three exphcit s,mple examples (for 
more examples see In particular we obtain a vanishing von Neumann entropy of a 

fermionic or a bosonic state containing the least possible amount of correlations. We believe 
that this settles an issue that has caused a lot of confusion regarding the use of von Neumann 
entropy as a n,easure of entanglement for .dentieal part.eles Q 



THE GNS CONSTRUCTION 

A general quantum system is usually described in terms of a Hilbert space T-L and linear 
operators acting thereon. Physical observables correspond to self-adjoint operators {O = 

: Ti ^ Ti). The probabilistic character of the theory is based on the notion of state, from 
which probabilities and expectation values can be computed. Generically, a state is described 
in terms of a density matrix p : H — ?■ "H, a linear map satisfying Tr p = 1 (normalization), 

= p (self-adjointness) and p > (positivity). For pure states, the additional condition 
p^ = p is required, so that p is of the form \ip){il)\ for some normalized vector \ip) G "H. 

Since the expectation value of an observable O is defined by {O) p = Tr(p(9), we can 
equivalently regard p as a linear functional O i— )■ {0)p from the space of operators to 
C. Moreover, since the space of all (bounded) operators on % forms an algebra €{1-1), it is 
possible to give a formulation of quantum physics which does not a priori make use of Hilbert 
spaces. Such a formulation was initially envisaged by von Neumann. The formulation due to 
Gel'fand and Naimark and further developed by Segal (GNS construction) led to the notion 
of an "abstract algebra of physical observables" , or C*-algebra. This construction (explained 



3elow) has played a very important role in quantum field theory 10| and statistical mechanics 
1^ . We propose to show that this approach is also very well-suited to deal with the problem 
described in the introduction. 
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We thus consider an abstract algebra A (playing the role of CiTi) above) that represents 
the physical observables. Since these observables are (not yet) acting on any Hilbert space, 
an abstract notion for the adjoint of an operator is required. This is provided by an operation 
("involution") a i— )■ a*. The algebra is assumed to contain an identity 1_4 and to be closed 
under products, linear combinations and under the involution. In this context, a state is 
defined as a linear functional u : A ^ C Again, since there is no Hilbert space, no density 
matrix appears at this stage. But from the interpretation of uj{a) = {a)^^ as the expectation 



value of a, the conditions of normalization uj{1X) = 1, reality uj{a*) = uj{a) and positivity 
u{a*a) > (for any a G A) are physically motivated properties that any state u must, by 
definition, satisfy. 

Given a quantum system defined by an algebra A and a state u, how do we recover the 
usual Hilbert space on which the algebra elements act as linear operators? Since A is an 
algebra, it is in particular a vector space, denoted here as A. Elements a E A regarded as 
elements of the vector space A are written as \a). Then, (3 & A will act on \a) G ^ as a 
linear operator by [5\a) := \[5a). A similar construction occurs when we study the regular 
representation of a group through its action on its group algebra [l^. 

In order for the vector space A to become a Hilbert space, an inner product is required. If 
we set = aj(a*/3), we obtain almost all properties of an inner product. In fact, reality 



and positivity can be used to show that (/3|a;) = (a;|/3) and also that {a\a) > 0. But it can 
happen that {a\a) = for some non-zero elements a. That is, there could be a null space 
AC of zero norm vectors: Af^ = {\a) G A\uj{a*a) = 0}. The solution to this problem is 
obtained by considering the quotient vector space A/Af^. Its elements are equivalence classes 
I [a]), with I [a]) equivalent to precisely when a — /3 G AC- In particular, if a G A^, then 
I [a]) = 0. The space A/M^^ has now a well-defined scalar product given by 

{[aWm = uj{a* P\ (1) 

independently of the choice of a from [a] and with no non-zero vectors of zero norm. Its 
closure is the GNS Hilbert space l-Lu)- In this way, one obtains a representation tt^^ of A on 



by linear operators 
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7r^(a)|[/3]) = IM). 

Partial Trace as Restriction - Consider a bipartite system 1-i = T-La ® "Hb, with a 
density matrix p. The description of Ha as a subsystem involves the reduced density matrix 
Pa, obtained through partial tracing over B. Using the language of algebras and states. 



we observe that the algebra corresponding to the joint system AB is given by ^ = C{l-L). 
Expectation values are computed using the state Up induced by p, (C)p = ^jJpiO) = Tr-^(p O). 
Corresponding to subsystem A, we can consider the subalgebra of "local" operators of 
the form (g) 1^, for an observable on T-La. We can then define a state Up^ : ~^ C 
which is the restriction Wpl^^ of ujp to defined by Up^Q^a) = Upi^a) if a G Aq. 

Now we observe that the reduced density matrix p^, obtained by partial tracing, gives 
rise to a state on subsystem A that is precisely the restriction of Up to 

cjp^o{K^lB) = TinAPAK). (2) 

Hence, partial trace and restriction give the same answer in this case. The importance of 
this observation lies in the fact that when "H is not of the form of a 'simple tensor product', 
partial trace is not a suitable operation. In contrast, if the system is described in terms 
of a state Up on an algebra A, it is still sensible to describe a subsystem in terms of a 
corresponding subalgebra Aq and of the restriction Up^ of Up to Aq. The GNS theory is 
well-suited for the study of Up^ for general algebras ^ ^■ 

von Neumann Entropy - The representation vr^ is in general reducible. This means 
that "Hoj can be decomposed into a direct sum of irreducible spaces: "H^ = 0j'Hi, where 
T!'uj{a)'Hi C Tii for all a E A. Let Pi : — t- Tii be the corresponding orthogonal projectors. 
These projectors can be used to construct a density matrix on the GNS space "H^ that 
yields the same expectation values as the original state u. The von Neumann entropy of p^j 
can then be evaluated in the standard way. The construction of p^^ goes as follows. 

First, we observe that the identity 1^ of A satisfies l^a = a for all a G as well 
as 1^ = 1^. This, together with ([T]), implies u{a) = Since the linear oper- 

ator TCuiia) is defined by 7r(^(a)|[/3]) = | [«/?]), we know that \[a]) = n (a) \[1 _/(]). It fol- 
lows that u{a) = ([l^]|7r,(a)|[l^]). Using = Ei.Pill^A]), = 
and from the orthogonality of the projectors, one obtains u{a) = Tr^^ {puTTuici)), where 
p^ = J2i^i\[^A]){[^A]\Pi- The von Neumann entropy of p^^ is then given by S{p^) = 
- l^i^og2 Pi, where pi = ||Pi . 

The crucial fact is that u is pure if and only if the representation tt^ is irreducible. In 
particular, the von Neumann entropy of u, S{u) = S{p^), is zero if and only if Ti^^ is 
irreducible. This property depends on both the algebra A and the state u. 

Consider now a subalgebra .4o C ^ of A. Let Uq denote the restriction to of a pure 
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state w on ^ jol- We can apply the GNS construction to the pair (^0,1^0) and use the von 
Neumann entropy of uq to study the entanglement emergent from restriction. 

Bipartite Entanglement from GNS - We now illustrate how to apply the GNS 
construction to entanglement. Consider "H = "Ha ®'Hb = (S) C^. The algebra of the full 
system is ^ = M2(C) (g) M2(C). Let us consider the normalized state vector (0 < A < 1): 
li'x) = V^|+, — ) + a/(1 — A)| — , +) , with corresponding state u on the algebra A: uj{0) = 
{HO\^x), OeA. 

Entanglement of I^^a) is to be understood in terms of correlations between "local" mea- 
surements performed separately on subsystems A and B. Measurements performed on A 
correspond to the restriction uq = u of u to the subalgebra Aq C A generated by el- 
ements of the form a ® I25 with a G M2(C). We obtain Wo(« ® I2) = (V^aIc^ ® 12|V^a) = 
A(+|q;|+) + (1 — A)(— |q;|— ). In accordance with ([2]), we have Ci;o(a® I2) = Tt^c^ {PaC(), where 
Pa = Tr^ \ipx){i^x\, namely, 

'■=(:.->) 

Now we perform the GNS construction based on the algebra ^0 — ^2(C) and the state 
uq. These are the data needed to describe subsystem A. For a G M2(C), we have coo{a) = 
Xaii + (1 — A)a;22- Now we consider ^0 as a vector space. This is just the assertion that 
M2(C) is, by itself, a vector space. From the explicit form of uq, one readily concludes that, 
as long as < A < 1, there are no null states. This means that the GNS space "H^g is just the 
four dimensional space of 2 x 2 matrices, endowed with the inner product {a\(3) = uola'^P). 
We can consider a basis of four 2x2 matrices defined as Cij = \i){j\ for i,j G {1,2}, 
where |1) = |+) and |2) = |— ). Then, for example, (eii|eii) = A and (e22|e22) = 1 — A. 
With an appropriate normalization and ordering of this basis, one checks that the operator 

( a o\ 

corresponding to a G ^0 is the 4x4 matrix 7r^o(a) = , showing in an explicit way 

y ay 

that the representation is reducible (the GNS-space splits as H^^^ = © C^). Following the 
prescription described above one obtains, for the density matrix, p^^^ = diag{A, 0, 0, 1 — A}. 
The identity ix'o(a) = Tr^^ o(p(^o7r^o(a)) is readily checked. 

From the explicit form of p^jQ we conclude that the entropy of the restricted state is 
S{ijJo) = —A log2 A — (1 — A) log2(l — A). This is precisely the entropy of the reduced density 
matrix p^ obtained by partial tracing. Recalling that a (pure) state of the full system is 



entangled with respect to a bipartition into subsystems if and only if S{pa) > 0, we have thus 
verified that our method reproduces the standard results in the case of bipartite systems. 
This is in fact a general result: 

For bipartite systems of the form Ha ®'Hb (pure case), the GNS construction yields a 
vanishing entropy for the restricted state precisely when the original state of the full-system 
is separable. Moreover, in the case of entangled states, the entropy computed via the GNS 
construction coincides with the von Neumann entropy of the reduced density matrix computed 
via partial trace and can therefore he used as an entanglement measure. 

We remark that, in the pure case, entanglement can also be characterized by the impos- 
sibility of writing the state as a product state. That is, if u is of the form u{a ® (5) = 
u}a{ol)^b{I3) for a {(3) any observable on subsystem A (B) and cja, ojb states on the corre- 
sponding subsystems, then a; is a product, or separable state, and it is not entangled. This 
observation makes clear that entanglement for mixed states can also be studied using our 
approach: If a mixed state ui,^ can be written as a convex combination of product states, 
then it is called separable, otherwise it is called entangled. 

SYSTEMS OF IDENTICAL PARTICLES 

In the case of identical particles, the Hilbert space of the system is not anymore of the 
tensor product form. Therefore, the treatment of subsystems using partial trace becomes 
problematic. In contrast, in our approach all that is needed to describe a subsystem is the 
specification of a subalgebra corresponding to the subsystem. Then, the restriction of the 
original state to the subalgebra provides a physically motivated generalization of the concept 
of partial trace, the latter not being sensible anymore. Applying the GNS construction to 
the restricted state, we can study the entropy emerging from the restriction and use it as a 
generalized measure of entanglement. 

Let = be the Hilbert space of a one-particle system. The fc-particle Hilbert 

space for bosons (fermions) is the symmetrized (antisymmetrized) /c-fold tensor prod- 
uct of 1-L^^\ To any one-particle observable A^^^ on 'H'^^\ we can associate the operator 

^(fe) (^(1) <g) 1^ . . . (g) 1^) + (1^ (g) ^(1) (g, . . . 1^) + . . . + (g, . . . (g, 1^ ^(1)) on n^^\ 

The operator A^^'> preserves the symmetries of 'HS^\ The map A^^"* — )■ A^^'> allows us to 
study subalgebras of one-particle observables. These constructions are most conveniently 

7 



expressed in terms of a coproduct A 15|. In fact, an approach based on Hopf algebras 15 1 
has the advantage that para- and braid-statistics can be automatically included. In what 
follows we use the simple coproduct A{g) = g ® g, g & U{d), linearly extended to all of 
CU{d). It gives the formula for A^''^ at the Lie algebra level. Physically, the existence of 
such a coproduct is very important. It allows us to homomorphically represent one-particle 
observables in the fc-particle sector. In the examples considered below, observables on such 
identical-particle systems can also be described in terms of creation/annihilation operators. 

In the following examples we will concentrate, for the sake of clarity, on systems of 
two fermions and two bosons (more examples will be presented in a forthcoming paper). 
However, our methods can be easily generalized to study many-particle entanglement. 

Two Fermions - Consider, as in [4i|, a one-particle space describing fermions with two 
external degrees of freedom (e.g. 'left' and 'right') and two internal degrees of freedom (e.g. 
'spin 1/2'). They are described by fermionic creation/annihilation operators al/\ with 
a standing for 'left', b for 'right' and A = 1, 2 for spin up and down, respectively. The single- 
particle space is therefore Ti^^^ = C^. The two-fermion space is given by "H*-^^ = /\^ C'^ C 
H^^^ 'H^'^^ (A denoting anti-symmetrization). T-L^'^'^ is generated from the "vacuum" 1^2) 
using pairs of creation operators. An orthonormal basis is given by the vectors ajaglfi), 
fol^gl^) si-iid a^^b\,\Q), with A, A' G {1,2}. The two-particle algebra A of observables is thus 
isomorphic to the matrix algebra Mq{(C). 

For \-ipg) = {cos9a[bl + sm9alb\)\ri) , the corresponding state cog is given by uje{a) = 
{iljg\a\ipg) for a E A. We now choose the subalgebra to be given by the one-particle 
observables corresponding to measurements at the left location. It is generated by 1_4, = 
a\aiala2, Na = a\ai-\-ala2 and Tj=i^2,3 = (1/2) CL\{cri)^''^' dy . Now we consider the restriction 
of ujg to ^0 and study the GNS representation corresponding to this choice. For < 9 < 11/2, 
the null space turns out to be spanned by |?2i2) and — A^^). Therefore, the GNS- 
space Tig is four-dimensional and spanned by and {| [Tj])}j=i_2,3- One may show that 

Tig = Til® ^2, with Til spanned by |[Ti -|- 1X2]) = \ [a{a2]) and |[a2a2]), and 7^2 spanned by 
|[a|ai]) and \ [Ti — iT2]) = \ [alai]). The two representations are isomorphic. Moreover, from 
the decomposition = |[ct2'^2]) + lictiCti]) of {[^a]) into these irreducible subspaces, we 

obtain the entropy S{9) = — cos^ 9 log2 cos^ 9 — sin^ 9 logg sin^ 9. 

For 9 = 0, the null space is spanned by 1^12), |1^ — a\ai), |a2a2), |oifl2)- The GNS-space 
T-Lq is and isomorphic to the above 7^2- Similarly, for = 7r/2 we find that the GNS-space 
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is isomorphic to the above 1-Li. Both GNS-spaces are irreducible, so that the corresponding 
Ci;o,o and ujil^ are pure states with zero entropy. 

This result should he contrasted with the entropy S = log92 obtained via partial trace for 
states with Slater rank one such as ojq^q, co^fi above (cf. l^, l2/ and references therein), that 
correspond to simple Slater determinants and, therefore, should not be regarded as entangled 
states. 

Two Bosons - Consider the one-particle space = with an orthonormal basis 
|e2), lea)}. The two-boson space "H^^^ is the space of symmetrized vectors in 'H^^^i^'H^^K 
It corresponds to the decomposition 3(8)3 = 6©3of SU{3). An orthonormal basis for 7^^^^ 
is given by vectors {|ej V ej)}jjg{i^2.3} where V denotes symmetrization (and the vectors are 
normalized). The two-boson algebra of observables A^"^^ is thus isomorphic to Mq{C). 

For the particular choice \4'e,<j>) = sin6'cos0|ei V 62) + sin6'sin0|ei V 63) -|- cos 6*163 V 63), 
the corresponding state is cjg,,/, defined by U0^^{a) = {ipe,(j>\o:\ipe,(j)) for any a E A. For the 
sake of concreteness, we choose ^0 to be given by those one-particle observables pertaining 
only to the one-particle states |ei) and |e2). 

We consider the restriction wg^^ l^^. The 6 representation under the SU{2) action on |ei) 
and 1 62), splits as 6 = 3©2©1. The subalgebra ^0 is given by block-diagonal matrices. Each 
block corresponds to one of the irreducible components in the decomposition 6 = 3 © 2 © 1. 
The dimension of ^0 is therefore 3^-1-2^ + 1^ = 14. 

The construction of the corresponding GNS-representation follows the same procedure as 
in the previous example. The von Neumann entropy as a function of the parameters is given 
by S{9,(f>) = — sin^ 6'[cos^ 01og2(sin 6* cos0)^ -|- sin^ 01og2(sin 6* sin0)^] — cos^ 6* log2(cos 6*)^. 



CONCLUSIONS 

The strong point of our approach is that it provides a precise, universal, and mathe- 
matically natural way to characterize and quantify entanglement for systems of identical 
particles. For many years it has been known that the von Neumann entropy based on par- 
tial tracing does not give the physically correct answer when applied to systems of identical 
particles. Different (i.e. non-universal) criteria have been developed which strongly depend 
on the statistics of the particles. In contrast, our approach is conceptually clear and ap- 
plies equally to any quantum system. It thus promises to resolve the controversy regarding 
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entanglement of identical particles [7|. 
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